function [x sign means] = bactrian(step, distance, sigma1, sigma2)
%BACTRIAN sum of 2 gaussians where the distance between gaussian's means.
%Gaussians are placed symmetrical to 0.
%   INPUT:
%       step - distance between sequential argument values
%       distance - distance between gaussians
%       sigma1 - standart deviation of first gaussian
%       sigma2 - standart deviation of second gaussian
%   OUTPUT:
%       x - argument values
%       sign - signal values
%       means - generated gaussians' means argument values
%       fmeans - generated gaussians' means function values 

x = -20:step:20;
% array of distribution means
means = [-distance; distance]/2;

sign = rightgaussmf(x, sigma1, means(1));
sign = sign + rightgaussmf(x, sigma2, means(2));